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The phase diagram of the quasi two-dimensional antiferromagnet BaNi2V20s is studied by specific 
heat, thermal expansion, magnetostriction, and magnetization for magnetic fields applied perpendic- 
ular to c. At fioH* c± 1.5 T, a crossover to a high-field state, where Tn(H) increases linearly, arises 
from a competition of intrinsic and field-induced in-plane anisotropics. The pressure dependences 
of T/v and if* are interpreted using the picture of a pressure-induced in-plane anisotropy. Even at 
zero field and ambient pressure, in-plane anisotropy cannot be neglected, which implies deviations 
from pure Berezinskii-Kosterlitz-Thouless behavior. 
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The study of quasi two-dimensional (2D) magnetic sys- 
tems [1] continues to be a focus of theoretical and ex- 
perimental investigations, motivated in large part by the 
discovery of high-temperature superconductivity in the 
quasi 2D cuprates. Further, the search for a magnetic 
system exhibiting true Berezinskii-Kosterlitz-Thouless 
(BKT) behavior, initially proposed for 2D XY magnetic 
systems [2j, has been elusive and has only been seen 
in superfluid and superconducting films [3]. Theoretical 
studies indicate that BKT behavior can also be expected 
for 2D Heisenberg systems with a small easy-plane XY 
anisotropy [4]. Two recent experimental papers suggest 
that BaNi2V2 08 may in fact be a physical realization 
of such a system [H|, 0]. BaNi2V208 has a rhombohe- 
dral structure (space group R3) and its magnetic proper- 
ties arise from a honeycomb-layered arrangement of spins 
S = 1 at the Ni 2+ sites. The quasi 2D properties are due 
to a strong antiferromagnetic Heisenberg superexchange 
J in the NiO honeycomb layers. Long range antiferro- 
magnetic ordering, which would be precluded in a purely 
2D Heisenberg system, sets in below the Neel tempera- 
ture T/v — 50 K [5] because of small additional energy 
scales, which we include in the following Hamiltonian: 
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The planar XY anisotropy Dxy — 1 meV is a factor of 10 
smaller than J [7] and confines the spins to lie within the 
honeycomb layers (easy plane). If this were the only ad- 
ditional term in Eq. (pQ), a true BKT transition could be 
expected within each 2D layer [4]. However, a real crystal 
is always three-dimensional (3D) and a very small inter- 
layer exchange J' ultimately leads to a crossover from 2D 
to 3D correlations, and then to a 3D ordering transition 
The value of J' is unknown for the present case 0; 



however, the extremely small signal at T/v in the specific 
heat [5] suggests that J' / J is very small (typically, J' / J 
is in the range 10 -2 -10 -6 in quasi 2D systems |l|). The 
in-plane anisotropy Dip is estimated by 4*10 ~ 3 meV Q 
and acts to align the spins along one of the three equiva- 
lent hexagonal easy a-axes [5|. The last term in Eq. (pQ) 
includes the effect of a magnetic field which also acts 
as an effective anisotropy y], @] . 

In this Letter, we study the (T, H) phase diagram of 
BaNi2V20s for magnetic fields applied within the honey- 
comb planes. The combination of specific heat, thermal 
expansion, magnetostriction, and magnetization allows 
us to calculate the uniaxial pressure dependences of T/v, 
J, and of a crossover field, iJ*, related to spin alignment. 
In particular, we find a strong anisotropy of the pressure 
dependences of T/v in the easy plane, which is related to 
the spin direction. This directly underlines the impor- 
tance of the in-plane anisotropy term Dip in Eq. (pQ) 
for establishing 3D long range order in BaNi2V2 08- Our 
analysis and comparison with literature indicates that 
the obtained (T, H) phase diagram may be generic to 
quasi 2D antiferromagnets with a degenerate number of 
easy axes and with J ^> Dxy ^> Dip, J' . It would be 
useful if the theoretical models for quasi-2D Heisenberg 
systems, which already consider XY anisotropy and inter- 
plane coupling [4], would also incorporate this in-plane 
anisotropy. 

Single crystals of BaNi 2 V 2 8 (42.8 mg) and 
BaNio.8Mgi.2V2 08 (4.85 mg) were grown in fluxes com- 
posed of BaC03, NiO, MgO, and V2O5, using AI2O3 cru- 
cibles. Specific heat was measured using a PPMS from 
Quantum Design Thermal expansion and magne- 
tostriction were measured using a high-resolution capac- 
itive dilatometer with temperature and field sweep rates 
of 20 mK/s and 0.5 T/min, respectively. The cell is ro- 
t at able, allowing both longitudinal and transverse mea- 
surements. Magnetization was measured using a MPMS 
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from Quantum Design. For all measurements, the mag- 
netic field was applied parallel to the easy plane. 

Fig. Q](a) shows the 0-T specific heat C p of BaNi 2 V2 08 
in a Cp/T versus T plot. An estimation of the phonon 
contribution C^ h is made using the specific heat of 
BaNio.8Mgi.2V2 08 pl|. The resulting magnetic specific 
heat C™ ag = C p - Cf is shown in the insert of Fig. [Q 
(a), where C™ ag /T is plotted as a function of T. The in- 
tegrated entropy ASmag — 24 J/mol K is roughly equal 
to the value 2i?ln3 ~ 18 J/mol K expected for the S = 1 
Ni 2+ ions. In this plot, there is little signature of the 
transition at T/v; rather, C^ 710,9 /T has a broad maximum 



at Tyy 



70 K which is attributed to the build-up of 



2D correlations and corresponds to the highest magnetic 
energy scale J. Fig. [T] (b) shows a blow-up of the C p 
anomaly related to the ordering at T/v- At H = 0, we 
clearly observe a tiny peak yielding T/v = 47.4 ±0.1 K, 
defined as the locus of the minimum of d(C/T)/dT [Tlj . 
The magnetic entropy ASn — 13 mJ/mol K gained at 
T/v is only a tiny fraction of the total magnetic entropy, 
the very small ratio ASn / ASm a g ^7* 10 -4 [ll[ being a 
consequence of the strong 2D Heisenberg character of the 
system |l2j. Application of a magnetic field in the easy 
plane leads, at 10 T, to the increases of T/v by about 3 
K and of ASn by a factor 4. 

The linear thermal expansivity a — (l/L)dL/dT is 
shown in Fig. [2](a-c) for T close to T/v and < fioH < 10 
T. Data for three configurations are presented: a longitu- 
dinal one, ai ( L || H || a), and two transverse ones, 
(L || a and H 1 (L,c)), and a c (L || c and H || a). At 
H = 0, T/v is characterized by the jumps Aa a > (Fig. 
[2](a-b)) and Aa c < (Fig. [21(c)). Application of a mag- 
netic field in the easy plane induces an increase of T/v, 
defined at the extremum of da/dT, in agreement with 
the specific heat data. The field induces a sign change 
of Aa for the longitudinal configuration (Fig. [2] (a)), 
while the sign of Aa does not change for the transverse 
configurations (Fig. [2](b-c)). 

Fig. [3] shows the values of T/v extracted from specific 
heat and thermal expansion. The key observation, which 
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FIG. 1: (color online) (a) Specific heat of BaNi2V20s and 
BaNio.8Mgi.2V208, plotted as C p /T versus T; the estimated 
magnetic specific heat of BaNi2V20s is plotted in the insert, 
(b) Cp/T versus T of BaNi 2 V 2 8 with H _L c. 




FIG. 2: (color online) (a-c) a— and c— axes thermal expansiv- 
ity of BaNi2V208 versus temperature, the field being applied 
parallel to the easy-plane, i.e. H _L c. (d-e) Schematics of the 
magnetic structure of the Ni 2+ ions in the plane and of the 
field-induced lattice distortion (below Tat). 



is independent of the direction of H within the hexagonal 
plane, is that Tn(H) is almost constant for /hqH < 2 T 
and increases linearly for ji^H > 2 T. The derivatives 
of isothermal magnetization and magnetostriction (Fig. 
[H (d)) both show broad peaks, without any measurable 
hysteresis, at n^H* ~ 1.5 T. As seen in Fig. El fioH* 
is nearly temperature independent, varying between 1.5 
T at 5 K and 2 T near T/v. We interpret H* as a spin- 
flop-like crossover field, above which the spins are aligned 
nearly perpendicular to the field (Fig. [2] (e)). Contrary 
to a first-order spin-flop transition, which occurs in a 
single-easy-axis system when H is applied parallel to the 
spins, a crossover is obtained because of the presence, at 
zero field, of three kinds of domains, in which the spins 
are orientated along equivalent hexagonal directions 0|. 
In addition to the collective rotation of the spins induced 
by the field, a domain alignment, by domain wall motion, 
probably plays a crucial role in the alignment of the spins. 
This crossover may additionally be broadened because of 
a tiny domain wall energy due to the quasi 2D nature of 
the magnetic exchange |l|. 

For H _L c, the alignment of the spins, which is con- 
trolled by the minimization of |H • S|, can be represented 
by an effective field-induced Ising in-plane anisotropy 
D e / p f (H.) 0, in which the easy-axis is perpendicular 
to (H,c). In the high-field AF' state (Fig. E]), i.e. for 
/loH > 2T, D e / p f (H) dominates over the intrinsic Dip 
term and the spins are almost perpendicular to H. The 
linear increase of Tn(H) observed in this regime is related 
to a reduction of the spin fluctuations along the direction 
of H, due to the field-induced anisotropy [3]. The high- 
field line extrapolates to H = at T/v,o = 46.6 ±0.1 
K, which is smaller than T N (H = 0) by A = 0.8 ± 0.2 
K. Theoretical support is needed to quantitatively relate 
T/v,o and the slope of Tjy(H) to the characteristic en- 
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FIG. 3: (color online) Phase diagram obtained with H J_ c. 
PM, AF, and AF' denote the paramagnetic, the low- and the 
high-field antiferromagnetic phases, respectively. 



ergy scales of the problem. Nevertheless, we speculate 
that T/v o is the ordering temperature of the system in 
the limit of no in-plane anisotropy, i.e. with Dip — > 
and H* — > 0, so that T/v,o is only controlled by J, J', 
and Dxy- In this picture, the hexagonal Ising-like in- 
plane anisotropy Dip stabilizes the long-range magnetic 
ordering, shifting Tjv,o upwards to T/v by A. Since the 
3D character of the long-range ordering was shown by 
neutron scattering [5], we interpret our phase diagram 
as resulting from a field-induced crossover from a 3D 
Ising-like long-range ordering controlled by the hexag- 
onal in-plane anisotropy D/p, to a 3D Ising long-range 
ordering controlled by an effective field-induced in-plane 
anisotropy D e /J (H). As long as there is a distribution 
of domains, i.e. for H < H*, T/v (if) is independent of H 
and is not controlled by D e /J (H). We speculate that our 
phase diagram - and its interpretation - may be generic 
to quasi 2D systems with more than one easy axes and 
with J > D X Y > D IP , J' pjjj|. 

In the following, we use our data, together with ther- 
modynamic relationships, to determine the uniaxial pres- 
sure dependences of T/v, J, and H*. The Ehrenfest rela- 
tionship dT^/dpi = AaiVTjy/ AC P , where Aa$ and AC P 
are the expansivity and specific heat jumps at the transi- 
tion, gives the uniaxial pressure dependences of T/v. Sim- 
ilarly, the uniaxial pressure dependences of J can be ob- 
tained by substituting T max ~ J for T/v in an Ehrenfest- 
type relation and considering the total magnetic signal 
for Aa and AC P . Appropriate scalings of specific heat 
and thermal expansion (at H = 0) are shown in Fig. [4](a- 
c) [16| , and the resulting pressure dependences are listed 
in TableU The hierarchy (l/J)dJ/dp < (l/T N )dT N /dp 
indicates that dTjy/dp is not controlled by J. We also 
note the field-induced sign change of dT^/dpa when 
H || a. In contrast, dTx/dp a and dT^/dpc remain pos- 
itive and negative, when H _L (p, c) and H || a, respec- 
tively, over the whole investigated H range. In Fig. [H 
(d), the derivative of magnetization dMfd(jiioH)^ and the 



magnetostriction coefficient a|J = (l/L a )dL a /d(fioH) 
measured in the longitudinal configuration, are shown 
at 5 K as a function of H. The data were scaled at the 
peak values to apply the generalized Ehrenfest relation 
dH*/dpi = VA\i/A(dM/d(n H)). The obtained p de- 
pendences of H * (Table [Tj) are much larger than the ones 
of T/v, which in turn are larger than those of J, and they 
are the largest for uniaxial pressures applied within the 
hexagonal plane. Since H* depends strongly on the zero 
field in-plane anisotropy, this suggests that the pressure 
dependence of the in-plane anisotropy controls the pres- 
sure dependence of H* . Further, in-plane pressure effects 
on H* have opposite signs parallel and perpendicular to 
the applied field. In the following discussion, we provide 
a simple explanation of the above correlations. We note 
that the c-axis pressure dependence of H* is much smaller 
than the in-plane pressure dependences of if* and may 
result from them via elastic coupling of the axes. 

In the high- field phase (AF' in Fig. [3]), the spins are 
aligned such that their direction is perpendicular to H. 
Associated with this spin-flop-like crossover is a macro- 
scopic distortion of the hexagonal plane, which develops 
below Tat and can be extracted from our dilatometry 
data (the integration of a(T) and A(T), in Fig. [2] (a-b)) 
and 2] (d), leads to the length changes). This in-plane 
distortion is such that the crystal contracts along the 
spin direction and expands perpendicularly to the spins 
(Fig. [2] (e)). Since presumably the magnetic domains 
are already distorted at low fields [17[ , the application of 
uniaxial pressure p in the easy-plane will tend to align 
the spins such that the contracted direction is parallel to 
p. Hence, the favored spin direction will be parallel to p. 
This effect is analogous to the magnetic field effect and 
can be described by adding an effective pressure-induced 
in-plane Ising anisotropy term D e /J (p) to the Hamilto- 
nian, the corresponding easy axis being parallel to p. 

The H and T evolutions of this in-plane distortion ulti- 
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FIG. 4: (color online) (a-c) Scalings of specific heat and ther- 
mal expansivity (with a minus sign for the c-axis) used to 
determine the zero-field uniaxial pressure dependences of Tn 
and J. (d) Scaling of dM / 'd(fioH) and Aa versus H at 5 K. 
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TABLE I: Normalized uniaxial pressure-dependences of J, 
Tjy, and H* obtained from our data using Ehrenfest relations. 
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mately govern the in-plane uniaxial p dependences of H* 
and Tat. p applied along a x , i.e. perpendicular to (H, c), 
will favor the alignment of the spins perpendicular to H, 
and thus will reduce H*. Using (l/H*)dH*/dp = -1.4 
kbar -1 , we estimate that 0.5 kbar along a x is enough to 
reduce H* to zero and to align the spins by stress alone. 
On the other hand, (l/H*)dH* /dp is positive for p ap- 
plied along a", i.e. parallel to H, which means that larger 
fields will be needed to align the spins. Magnetic ordering 
will be favored if both the field- and the pressure-induced 
anisotropics D e /J (H) and D e /J (p) act cooperatively to 
align the spins along the same axis [18[ . This explains the 
positive value of dTjy /dpj^ at 10 T. In contrast, D e /J (p) 
acts against D e /J (Yi) when p is applied along a", i.e. 
parallel to H, which leads to dT N /dpl < at 10 T. 

In conclusion, the (T, H) phase diagram of BaNi2V208 
was studied with H _L c. The obtained field- and uniaxial 
pressure-dependences of T/v were interpreted as the con- 
sequence of effective field- and pressure-induced in-plane 
anisotropics. Our data clearly demonstrate the impor- 
tance of the in-plane anisotropy for establishing 3D long- 
range order in BaN^X^Os- A search for BKT behavior 
should be restricted to a temperature region above T/v, 
where in-plane anisotropy no longer influences the growth 
of the 2D-XY correlations. In the limit of Djp = 0, 
we extrapolated a Neel temperature Tn,o = 46.6 K, 
which provides an upper estimate of the BKT temper- 
ature Tbkt [19]. Since the Ising-like anisotropy Dip 
leads to an increase of T/v by almost 1 K, 2D-XY BKT 
behavior may be valid for temperatures at least several 
degrees higher than T/v- These results are expected to 
be useful for future studies of BKT behavior. 
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